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Combining strong magnetic anisotropy with strong exchange interaction is a long standing goal 
in the design of quantum magnets. The lanthanide complexes, while exhibiting a very strong ionic 
anisotropy, usually display a weak exchange coupling, amounting to only a few wavenumbers. Re¬ 
cently, an isostructural series of mixed Ln^^-Nj^-Ln^"*" (Ln = Gd, Tb, Dy, Ho, Er) have been 
reported, in which the exchange splitting is estimated to reach hundreds wavenumbers. The micro¬ 
scopic mechanism governing the unusual exchange interaction in these compounds is revealed here 
by combining detailed modeling with density-functional theory and ab initio calculations. We find 
it to be basically kinetic and highly complex, involving non-negligible contributions up to seventh 
power of total angular momentum of each lanthanide site. The performed analysis also elucidates 
the origin of magnetization blocking in these compounds. Contrary to general expectations the 
latter is not always favored by strong exchange interaction. 


INTRODUCTION 

The effects of strong magnetic anisotropy, tradition¬ 
ally investigated in magnetic insulators, especially, in /- 
electron systems [1-3], recently attracted renewed inter¬ 
est in connection with molecular magnetic materials [4]. 
The investigation of molecular nanomagnets gave birth 
to new objects such as single-molecule magnets (SMMs) 
[5, 6 ] and single-chain magnets [7], and initiated studies 
in the domain of molecular spintronics [ 8 , 9] and quantum 
computation [10-12]. Among them, in the last years the 
accent moved towards lanthanide complexes which have 
already demonstrated several exciting properties [13-18]. 

The key feature of lanthanide ions in materials is their 
strong magnetic anisotropy caused by strong spin-orbit 
coupling effects [19], which often leads to highly axial 
ground and low-lying excited doublet states even in the 
lack of axial symmetry [20]. Due to small radius of elec¬ 
tronic /-shells, the exchange interaction in lanthanide 
complexes is much weaker than the crystal-field splitting 
on lanthanide ions [ 21 ]. As a result, only individual dou¬ 
blet states on lanthanide sites, described by pseudospins 
s = 1/2, participate in the magnetic interaction. The lat¬ 
ter is described by a Hamiltonian bilinear in pseudospins 
(si and S 2 ) in the case of two interacting lanthanide ions, 
or a pseudospin (si) and a true spin ( 5 ' 2 ) in the case of 
a lanthanide ion interacting with a transition metal or a 
radical when the spin-orbit coupling in the second site is 
negligible. For strongly axial doublet states on the lan¬ 
thanide sites (Ln) these Hamiltonians basically become 
of Ising type [ 22 ]: 

^hn—S — \J^1ziS2z2^ ( 1 ) 

either collinear {zi || Z 2 ) or non-collinear (zi j/ Z 2 ) 
depending on geometry [ 21 ] and details of interaction. 
[23] The exchange parameter is contributed by mag¬ 
netic dipolar and exchange interaction between the sites, 
J' = J7dip + Jexch, the former being usually stronger in 
net lanthanide complexes [ 21 ]. 


This paradigm was recently challenged by a series 
of N 2 ~-radical bridged dilanthanide complexes [K(18- 
crown- 6 )]{[(Me 3 Si) 2 N](THF)Ln} 2 (/r-??^:r 7 ^-N 2 ) (Ln = 
Gd ( 1 ), Tb ( 2 ), Dy ( 3 ), Ho ( 4 ), Er ( 5 ), THF = tetrahy- 
drofuran), shown in Fig. la [14, 24]. In some of these 
compounds the exchange interaction was found to be two 
orders of magnitude stronger than in any known lan¬ 
thanide system. This is of the same order of magnitude 
as the crystal-held splitting of J-multiplets on the lan¬ 
thanide sites, implying that the picture of exchange in¬ 
teraction involving individual crystal-held doublets, Eq. 
(1), is no longer valid for these compounds. Moreover, 
the terbium complex from this series exhibits a magnetic 
hysteresis at 14 K and a 100 s blocking time at 13.9 K 
(one of the highest blocking temperatures among exist¬ 
ing SMMs [24]), suggesting a possible implication of the 
giant exchange interaction in this SMM behavior. 

The purpose of the present work is to reveal the mech¬ 
anism of giant exchange interaction and the origin of the 
magnetization blocking of the series of the complexes 
based on adequate theoretical treatment. We apply 
an approach combining ab initio and density-functional 
theory (DET) calculations with microscopic model de¬ 
scription to unravel the nature of this exchange interac¬ 
tion. We also elucidate the origin of blocking barriers in 
these compounds and discuss the effect of strength of ex¬ 
change interaction on magnetization blocking in strongly 
anisotropic complexes. 

RESULTS 

Origin of giant exchange interaction 

To understand the origin of such strong exchange in¬ 
teraction, we consider the simplest complex of the series, 
the gadolinium one. In this system the isotropic spins 
of Gd^+ ions (S'cd = 7/2) interact with the radical spin 
of the bridge (^Na = 1/2) via Heisenberg exchange 
interaction, iJex = Z[j=i ,2 “2 JHeisSGd(*) • Sns: described 
by a single parameter JB_eis due to the inversion symme- 
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FIG. 1. Molecular structure of Tb complex 2 and magnetic susceptibility in the series 1-5. a, Colors’ legend for the balls: violet, 
Tb; blue, N; red, O; green. Si; grey, C. The hydrogen atoms are omitted for clarity. The violet dashed lines show the orientation 
of the main anisotropy axes of Tb ions in their ground doublet state, whereas the green dashed line shows the orientation of the 
main anisotropy axis in the ground exchange Kramers doublet. The violet arrows show the orientation of the local magnetic 
moments on Tb ions, and the blue arrow on the radical, in the ground exchange Kramers doublet, b, Experimental (symbols) 
and ab initio calculated (lines) temperature-dependent powder magnetic susceptibility (x) for 1-5. The experimental data were 
upscaled by 3, 3, 1 % for 2, 3, and 5, respectively, and were downscaled by 2 % for 4. The magnetic susceptibility curves 
were calculated following the way they have been measured [14, 24], as /H at H— 1 T, averaged over all directions of 

magnetic field H relative to molecular frame. For the computational methodology of the magnetic axes and xT, see Refs. [25] 
and [26], respectively. 


try of the complex (Fig. la). Broken-symmetry DFT 
calculations [27] give the value Jneis = —21.4 cm“^ in 
close agreement with the experimental one, TTneis = —27 
cm“^ [14], and the previous DFT calculations [28, 29]. 

To get insight into the mechanism responsible for the 
obtained huge value of Jneis, we projected a series of 
DFT calculations into the effective tight-binding and 
Hubbard models acting in the space of interacting mag¬ 
netic orbitals of two Gd ions and the radical (see the Sup¬ 
plemental Material for details). Because of the D 2 h sym¬ 
metry of the exchange core (Fig. 2a), the antibonding 
TT* orbital accommodating the unpaired electron of 
radical overlaps with only one of the 4/ orbitals on each 
Ln site, the xyz one (Fig. 2b). The corresponding trans¬ 
fer parameter t was derived for the Gd complex as t = 
1407 cm“^. The value of t is obtained large because the 
radical’s magnetic orbital tt* resides on nearest-neighbor 
atoms (nitrogens) to both lanthanides. Most important, 
this orbital is found to lie higher than the 4:fxyz orbitals 
by as much as A = 5.2 x 10"^ cm“^ (Fig. 2b). Because 
of this huge energy gap, small electron promotion energy 
is expected for the electron transfer from the tt* to the 
'ifxyz orbitals: the Goulomb repulsion energy between 
the transferred electron and the / electrons is cancelled 
at large extent by A. On the other hand, because of 
the same large gap A, the promotion energy of electron 
transfer from 4/ to tt* orbital is at least one order of mag¬ 
nitude larger. Therefore, the contribution of this process 
to the exchange coupling can be neglected. Indeed, our 


analysis using the Hubbard model gives the experimental 
Jueis for the Gd complex with (averaged) promotion en¬ 
ergy of G = 8872 cm“^, a value many times smaller than 
typical “Hubbard D” in metal complexes [30]. Taking 
into account only the dominant virtual electron transfer, 
(4/)7(7r*)i ^ (4/)®(7r*)0 ^ {4fy{TT*)\ the kinetic con¬ 
tribution to the Gd^^-N®” exchange parameter is written 
in a good approximation as —2t‘^/U [31, 32]. 

Compared to this mechanism, the other contributions 
such as the direct exchange, the delocalization of un¬ 
paired electron of into the empty 5d orbitals of Gd®'*' 
(Goodenough’s mechanism [33]), the spin polarization 
and the magnetic dipolar interaction between Gd®“'" ions 
are expected to be 1 - 2 orders of magnitude smaller. The 
reason is that all these contributions are expected to be 
of the same order of magnitude as in other lanthanide- 
radical compounds. Indeed, the direct exchange integral 
depends only on the shape of the 4/ an radical’s orbitals, 
which is not expected to be much different from other 
complexes. The Goodenough’s contribution arises from 
higher (third) order of the perturbation theory compared 
to the usual kinetic exchange, and involves the excita¬ 
tion energy into a higher 5d orbital on the Ln site. Both 
these contributions are usually neglected unless the con¬ 
ventional kinetic exchange appears to be small [31, 32]. 
The spin polarization mechanism starts to play a role 
when the ligand bridging the magnetic centers contains 
a spectrum of low-lying orbital excitations, which is cer¬ 
tainly not the case of . As for magnetic dipolar inter- 
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action, it is estimated for Gd^^-N^” to be ~ 0.25 cm“^. 

The same physical situation is realized in the other 
complexes of the series. As Table I shows, the transfer 
parameters only slightly decrease with the increase of Ln 
atomic number. On the other hand, the gap A between 
the 4/ and the tt* orbital levels is obtained as huge as in 
the Gd complex (Table I), leading again to small promo¬ 
tion energy and, consequently, to the dominant role of the 
kinetic mechanism in the Ln^“'"-N 2 ~ exchange coupling of 
complexes 2 - 5 . Given the small change of t through 1 - 5 , 
the strong variation of the strength of exchange interac¬ 
tion in this series of complexes, testified by the experi¬ 
mental magnetic susceptibilities (Fig. lb), is expected to 
be due to the variation of the promotion energy. 


Anisotropic exchange interaction 


Contrary to the Gd complex, the other members of the 
series are characterized by strong magnetic anisotropy on 
the Ln sites induced by the crystal-field (CF) splitting of 
their atomic J multiplets. These CF-split multiplets are 
described by multi-configurational wave-functions, there¬ 
fore, they should be treated by explicitly correlated ab 
initio approaches [25, 26] rather than DFT. The ab initio 
fragment calculations show that the CF split J multiplet 
on Tb^+ ion (« 700 cm“^) is of the same order of mag¬ 
nitude as the estimated isotropic exchange splitting in 
1 (« 400 cm“^). Therefore, in sharp contrast with the 
common situation in lanthanides, the exchange coupling 
in the anisotropic 2-5 does not reduce to the interaction 
between individual (lowest) CF doublets on Ln sites with 
the S' = 1/2 spin of the radical, Eq. (1), but will intermix 
the entire CF spectrum arising from the ground atomic J 
multiplet at lanthanide ions. Then, such exchange inter¬ 
action should be formulated in terms of the total angular 
momenta (i = 1, 2) on the lanthanide sites. 

Extending the Anderson’s superexchange theory [31, 
32] to strong spin-orbit coupled systems, the tensorial 
form of the kinetic (covalent) interaction has been re¬ 
cently derived from the microscopic electronic Hamilto¬ 
nian [34]. The kinetic interaction between the lanthanide 
and radical centers contains besides the exchange part 
(Hex) also the Ln^^-N^” covalent contribution (arising 
from Ln^+-N 2 ~ electron delocalization) to the CF split¬ 
ting at the Ln^+ sites 


^cf = X! X! '^'=900 


i—1,2 kq 


-^ex — ^ ^ ^ ^ J^kqlq' 
i—1,2 kqq' 


oUh)i 
oliJ) ’ 

ol{j)s ■ 


( 2 ) 

( 3 ) 


Here, / and Sq are the unit and the spin operators, re¬ 
spectively, of the radical’s spin S = 1/2, 0^(J) are the 
Stevens operators [35] of rank k and component q, and 
fTfegoo and Jkqiq' are the exchange parameters [34]. The 
Stevens operator 0^(J) is a polynomial of Jq, (a = x, y, z) 


of /cth degree, in which [ql (= zLq) corresponds to the or¬ 
der of J± (= Jx ± iJy)- The maximal rank of fc is 7 
for the considered Ln^’*' ions, whereas the maximal [gl 
is 5 in the present case because only the '^fxyz magnetic 
orbitals at the lanthanide sites contribute to the kinetic 
exchange. The summation over k in Eqs. (2) and (3) 
is confined to even and odd ranks, respectively, which 
is required by the invariance of these Hamiltonians with 
respect to time-inversion. As it is seen from the form of 
these Hamiltonians, H/f only contributes to the CF split¬ 
ting of J multiplets on individual metal sites, whereas 
Hex describes the interaction between powers of total 
angular momenta at the metal sites with components of 
spin S = 1/2 of the radical. For comparison, the 
weak anisotropic exchange interaction between two spins 
(pseudospins) is described by the exchange Hamiltonian 
Si •D-S 2 , where D is the 3 X 3 exchange matrix, contain¬ 
ing one isotropic, Hve symmetric anisotropic and three 
antisymmetric (Dzyaloshinsky-Moriya) exchange param¬ 
eters [36]. This Hamiltonian corresponds to the first 
rank contribution (fc = 1) in Eq. (3), where Jiqiq' 
are just the nine components of the above exchange 
matrix D. The expression for the exchange parame¬ 
ter Jkqk'q' [34] includes all virtual electron transfer pro¬ 
cesses, (4/)”(7r*)^ —> (4/)”+^(7r*)° —> (4/)”(7r*)^, where 
n is the number of 4/ electrons in Ln^+. The multiplet 
electronic structure of Ln^“'" is fully included in the elec¬ 
tron promotion energy C/q -I- AH^ and the wave functions 
of the intermediate states, where by Hq we further denote 
the smallest promotion energy, a numbers the interme¬ 
diate J-multiplets, and AH^ is the excitation energy of 
the multiplet a with respect to the ground one in Ln^’*'. 

The highly complex tensorial form of the exchange 
Hamiltonian is inevitable for orbitally degenerate sys¬ 
tems with strong spin-orbit coupling, as was pointed out 
long time ago [37, 38]. Although all exchange parameters 
•Jkqk'q' are in principle required for adequate description 
of the exchange interaction, it is hardly possible to ex¬ 
tract a sufficient large number of them from experiment 
in a unique way. However, once Jkqk'q' are expressed 
via microscopic electronic parameters [34], the latter can 
be determined from up-to-date quantum chemistry cal¬ 
culations. Thus the transfer parameter t is obtained here 
from DFT calculations, expected to be accurate enough 
[39, 40], whereas the excitation energies AHq, and the CF 
states are obtained by fragment state-of-the-art ab initio 
calculations including spin-orbit coupling [25, 26]. The 
only parameter that might be inaccurate when extracted 
from DFT or ab initio calculations is Uq. Indeed, the 
former gives at most an averaged value over multiplets 
U and the latter systematically overestimates it due to 
insufficient account of dynamical correlation. 

In this way we construct the full microscopic Hamilto¬ 
nian, H = Hcf-l-H/£-|-Hex, containing only one unknown 
parameter Uq, where Hcf is the ab initio CF Hamiltonian 
for mononuclear Ln fragments (see Supplemental Materi¬ 
als for details). Diagonalizing this Hamiltonian, the mag¬ 
netic susceptibility y for the entire series of compounds 
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FIG. 2. a, Exchange core Ln^"'"-N 2 -Ln^^ in the complex corresponding to D 2/1 symmetry, b, Magnetic orbitals in 1 obtained 
from DFT calculations. Only the / orbital involved in the kinetic exchange mechanism is shown. 


TABLE I. Transfer parameters t, energy gaps A between the 4/ and the tt* orbital levels, minimal electron promotion energies 
Uq (all in cm“^), g-factors and angles between the magnetic moments on Ln®'*' and {6) in the ground exchange KD, and 
blocking barriers Ebarrier (cm~^) for complexes 1-5. For Ebarrier, both the experimental (exp.) [14, 24] and present (calc.) data 
are shown. 



1 (Gd) 

2 (Tb) 

3 (Dy) 

4 (Ho) 

5 (Er) 

t 

1407 

1333 

1322 

1311 

1270 

A 

5.20 X lO'^ 

5.74 X 10'* 

5.80 X 10* 

5.73 X 10* 

5.78 X 10* 

Uo 

8500 

4600 

6500 

7400 

12200 

9x 

2.2 X 10“^ 

7.6 X 10"® 

2.2 X 10“® 

4.7 X 10"® 

1.3 X 10“® 

9y 

3.7 X 10“2 

1.1 X 10"® 

7.0 X 10“® 

1.2 X 10"^ 

1.6 X 10"® 

9^ 

25.6 

33.6 

37.5 

36.2 

32.1 

6 

0 

0 

0 

2.5° 

CO 

0 

2.6° 

6.2° 

-E'barrier (^^P-) 

- 

227 

123 

73 

36 

-E'barrier (calc.) 

- 

208 

121 

105 

28 


has been simulated as described elsewhere [26]. Figure 
lb shows that the experiment is well reproduced for the 
values of minimal promotion energy Uq listed in Table I. 
The calculated exchange parameters for the series of the 
complexes are shown in Table II. We can see from the ta¬ 
ble that the exchange interaction involves non-negligible 
contributions up to the rank k = 7. 

The low-lying exchange spectrum for the Tb complex 
is shown in Fig. 3a. The ground (1±) and the first 
two excited (2±,3±) exchange Kramers doublets (KDs) 
mainly originate from the ground CF doublets on the Tb 
ions (94 %, 87 %, and 88 %, respectively). However, the 
third and fourth excited exchange KDs (4±, 5±) repre¬ 
sent almost equal mixtures of the ground and the first 
excited CF doublets on the Tb^"*" sites. This is remark¬ 
able because the mixed CF states are separated by 166 
cm“^ (Fig. 3a). Similar scenario is realized in 3 and 4, 


whereas in 5 the exchange interaction and the resulting 
mixing of CF states is relatively weak. The magnetic 
structure of the ground exchange KD is shown in Fig. 
la. The magnetic moments on Tb^"*" sites are parallel 
due to inversion symmetry and almost coincide with the 
directions of the main magnetic axes in the ground local 
KDs (Fig. la). The magnetic moment of the radical, cor¬ 
responding to isotropic S = 1/2, is rotated with respect 
to the magnetic moments on Tb sites by small angle 9 
(Table I) due to the non-Heisenberg contributions to the 
exchange interaction [23]. 

One may notice that the dominant first rank term of 
the exchange interaction is of isotropic Heisenberg type 
despite the strong spin-orbit coupling in Ln^+ ions (Ta¬ 
ble H). This looks surprising because even weak spin- 
orbit coupling makes the first-rank exchange interaction 
anisotropic [36]. The analysis of the expression for the 
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FIG. 3. The low-lying exchange spectrum and the magnetization blocking barrier in 2. a, The violet bold lines show the 
CF levels on Tb ions, the green bold lines show the low-lying exchange levels. Each exchange level is placed according to the 
projection of its magnetic moment on the main magnetic axis of the ground exchange doublet (green dashed line in Fig. la). The 
exchange levels with the same number are two components of the corresponding KD. The thin dashed lines show the admixed 
CF states on Tb sites to the exchange states in percent (only admixtures > 10% are shown). The number accompanying the 
red line is the average magnetic moment matrix element (in hb) between the components of the lowest exchange KD; the rate 
of QTM in the ground exchange state is proportional to its square. The red arrows denote the relaxation path outlining the 
barrier of reversal of magnetization, with the same meaning of the corresponding numbers (see the text for more details), b. 
The magnetization blocking barrier for 2 calculated in the absence of the admixture of excited CF states on Tb sites to the 
ground one via the exchange interaction. 


first-rank exchange parameters J\qiq' [34] shows that 
they are in general of non-Heisenberg type, whereas the 
present case is the only possible exception (see Supple¬ 
mental Material). Indeed, the isotropy of the first-rank 
exchange contribution requires involvement of only / or¬ 
bitals with the projections m = ±2. This can only arise 
for high symmetry of the exchange bridge (Fig. 2a) and 
for situations with one single electron transfer path, as 
in the present case. If any other orbital (or more of 
them) contribute to the electron transfer, the first-rank 
exchange interaction becomes strongly anisotropic. 


Magnetization blocking barriers 

Table I shows that the transverse g-factors {gx and gy) 
in the ground exchange KD, the squares of which char¬ 
acterize the rate of quantum tunneling of magnetization 
(QTM) [21], are the largest for 4 and the smallest for 
2 and 3 complexes. This explains why large magnetiza¬ 
tion hysteresis is seen at low temperatures in the latter 
two compounds, while not seen at all in the former and 
only weakly observed in the complex 5 [14, 24]. The 
path characterizing the activated magnetic relaxation in 
high-temperature domain is shown for the Tb complex 
in Fig. 3a by blue arrows. The height of the activation 
barrier ifbarrier corresponds to the first excited exchange 
KD, because its two components (2± in Fig. 3a) are 
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TABLE II. Calculated exchange parameters Jhq\q> (cm for 
the complexes 1- 5. 


k 

q 

<?' 

1 (Gd) 

2 (Tb) 

•^kqlq' 

3(Dy) 

4 (Ho) 

5 (Er) 

1 

0 

0 

94.9 

95.8 

70.8 

55.4 

24.2 

1 

±1 

Tl 

-94.9 

-95.8 

-70.8 

-55.4 

-24.2 

3 

0 

0 

0.0 

13.4 

-10.6 

-4.4 

5.0 

3 

±1 

Tl 

0.0 

8.2 

-6.5 

-2.7 

3.0 

3 

±3 

±1 

0.0 

10.6 

-8.4 

-3.5 

3.9 

5 

0 

0 

0.0 

17.0 

-16.0 

-1.6 

4.2 

5 

±1 

Tl 

0.0 

-12.8 

8.4 

6.8 

-6.1 

5 

±3 

±1 

0.0 

-2.5 

7.5 

-7.6 

3.4 

5 

±4 

0 

0.0 

5.7 

-0.8 

-7.5 

5.0 

5 

±5 

=Fl 

0.0 

-13.5 

11.5 

3.2 

-4.4 

7 

0 

0 

0.0 

0.3 

-3.3 

4.6 

-2.3 

7 

±1 

Tl 

0.0 

-0.2 

2.5 

-3.5 

1.7 

7 

±3 

±1 

0.0 

0.2 

-2.2 

3.0 

-1.5 

7 

±4 

0 

0.0 

-0.4 

5.1 

-7.1 

3.5 

7 

±5 

Tl 

0.0 

0.6 

-7.2 

10.0 

-5.0 


connected by a large magnetic moment matrix element 
which causes a large temperature-assisted QTM. Block¬ 
ing barriers of similar structure (Fig. 3a) arise in 3 and 4 , 
their calculated activation energies comparing well with 
the experimental ones (Table I). 

The unusually large matrix elements between the 
ground and the first excited exchange KDs are entirely 
due to the exchange mixing of the ground and the first ex¬ 
cited CF doublets on the Ln sites. Indeed, if one quenches 
the exchange admixture of excited CF doublets to the 
ground ones, this matrix element becomes three orders 
of magnitude smaller (Fig. 3b). Then the activated 
relaxation will proceed via a higher exchange doublet, 
thereby doubling the height of the blocking barriers (Fig. 
3b). Thus in the case of very strong exchange interac¬ 
tion, which is able to intermix the CF states on Ln sites, 
the axiality of the ground and excited exchange doublets 
is diminished dramatically and the blocking barriers do 
not exceed the energy of the hrst excited exchange KD. 
In other words, the strength of exchange interaction after 
reaching a certain value starts playing a destructive role 
for the magnetization blocking. Therefore, to exploit the 
effect of strong exchange interaction for achieving high 
magnetization blocking, an even stronger axial CF field 
on the Ln sites, precluding the exchange admixture of 
excited CF states, seems to be indispensable. 

DISCUSSION 

The mixed lanthanide complexes 1-5 investigated in 
this work are unique because they show an exchange in¬ 
teraction up to two orders of magnitude stronger than in 
conventional lanthanide complexes. Due to such strong 


exchange interaction, a qualitatively new situation arises 
when the exchange coupling starts to intermix the CF 
multiplets on the Ln sites. In all previous lanthanide 
complexes only the ground CF doublets on Ln sites were 
involved, which led to conventional Ising-type exchange 
interactions. In the present case, due to the involvement 
of all CF doublets belonging to the atomic J-multiplet, 
the exchange interaction becomes highly complex, requir¬ 
ing a tensorial description and involving many parame¬ 
ters. 

By combining DFT and ab initio calculations with the 
microscopic modeling of the exchange interaction, we 
were able to unravel the mechanism of giant exchange 
interaction in these complexes. This exchange interac¬ 
tion is found to be kinetic and highly complex, involving 
non-negligible contributions up to seventh power of total 
angular momentum of each Ln site. Based on the calcu¬ 
lated exchange states, the mechanism of the magnetiza¬ 
tion blocking is revealed. Contrary to general expecta¬ 
tions the latter is not always favored by strong exchange 
interaction. The accuracy of our approach is proved by 
the close reproduction of experimental magnetic suscep¬ 
tibility and magnetization blocking barrier for all inves¬ 
tigated compounds. 

The theoretical analysis proposed in this work opens 
the way for the investigation of highly complex exchange 
interaction in materials with strongly anisotropic mag¬ 
netic sites. Given the large number of involved exchange 
parameters and the obvious difficulties of their experi¬ 
mental determination, such an approach can become a 
powerful tool for the study of magnetic materials of pri¬ 
mary interest. 


METHODS 

DFT calculations. 

All DFT calculations were carried out with ORCA 
3.0.0. program [41] using the B3LYP functional and 
SVP basis set. Scalar relativistic effects were taken into 
account within Douglas-Kroll-Hess Hamiltonian. The 
isotropic exchange parameter for the complex 1, fTneisj 
was derived by applying the broken-symmetry approach 
[27]. The obtained Jneis was divided by 2 to account 
for its overestimation due to the self-interaction error 
[42, 43]. The 4/ and the tt* orbital levels and the trans¬ 
fer parameters t for all complexes 1-5 were derived by 
projecting the Kohn-Sham orbitals onto a tight-binding 
model. The averaged promotion energy U for the com¬ 
plex 1 was derived by reproducing the energy difference 
between the high-spin and the broken-symmetry DFT 
states with a Hubbard model. 
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Ab initio calculations. 

Energies and wave functions of CF multiplets on Ln^+ 
sites in 1-5 have been obtained from fragment ab initio 
calculations including the spin-orbit coupling, using the 
quantum chemistry package Molcas 7.8 [44], The calcu¬ 
lations have been done for the experimental geometry of 
the complexes, in which one of the two Ln^“'" ions was 
replaced by an isovalent closed-shell La^'*' ion. The total 
number of electrons was reduced by unity in order to have 
a closed-shell electronic configuration on the dini¬ 
trogen bridge. To simulate the electrostatic crystal field 
from the removed radical’s electron, two point charges 
of —0.5e were added on the nitrogen atoms. For this 


structural model of a Ln fragment, the complete active 
space self-consistent field (CASSCF) approach was used 
including all seven 4/ orbitals of the Ln atom in the ac¬ 
tive space. The spin-orbit interaction was treated with 
the module SO-RASSI and the local magnetic properties 
were calculated with the SINGLE_ANISO module of Mol¬ 
cas [45]. Exchange energy spectrum and magnetic prop¬ 
erties of the investigated polynuclear compounds were 
calculated using the POLY_ANISO program [26, 45], 
modified to treat the general form of exchange interac¬ 
tion, Eqs. (2), (3), within the kinetic exchange mecha¬ 
nism. 

For further details, see Supplemental Material. 
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Supplemental Materials 
for 

“Giant exchange interaction in mixed 
lanthanides” 

This material contains: 

1) DFT based derivations of the 4/ and the tt* orbital 
levels and of the transfer parameters t for all complexes 

1 - 5 ; 

2) Fragments ab initio calculations of the energies and 
wave functions of CF multiplets on Ln^+ sites in 1 - 5 , 
and calculations of atomic multiplets of the correspond¬ 
ing Ln^+ ions; 

3) The calculation of the exchange spectra are described; 

4) The analysis of the first rank exchange parameters. 


I. DFT CALCULATIONS 


A. Extraction of the transfer parameter t for 1-5 


In order to derive the transfer parameters between 
the 4/ orbital and the tt* orbital of the bridging N2, 
the Kohn-Sham levels are projected into tight-binding 
Hamiltonian: 


<7 


■ 2 

.i=l 


t I 




4" 



where i(= 1,2) is the index for the Ln^+ site in the com¬ 
plex, site is described by the type of the magnetic 
orbital tt*, 7 is the orbital component xyz, a =t,i is 
the projection of spin operator, tf and e,r*(= £/ + ^) 
are one electron orbital levels of the 4/ orbital and the 
TT* orbital, respectively, t is the transfer parameter be¬ 
tween the 4/ and the tt* orbitals, (c) is an electron 
creation (annihilation) operator, and h is a number op¬ 
erator. The subscripts of the creation, annihilation, and 
number operators indicate the site, the orbital index for 
only lanthanide site, and spin projection. Because of the 
D 2 h symmetry of the magnetic core part, only one 4/ or¬ 
bital {^fxyz) overlaps with the tt* orbital (Fig. 2b in the 
main text). Therefore, we only include the ‘ifxyz orbital 
for each lanthanide site in the model Hamiltonian. 

Diagonalizing the tight-binding Hamiltonian (SI), the 
one-electron levels are obtained as 



(S2) 

e/.s = ^ - V+ 8f^) , 

(S3) 

e-TT* ~ 2 ^, 

(S4) 


orbital levels with the DFT calculations, we obtain pa¬ 
rameters £/, t, and A. 

The highest occupied Kohn-Sham orbital for the down 
spin in the low-symmetry DFT solutions correspond to 
the TT* orbital. On the other hand, 4/ atomic orbitals 
contribute to many Kohn-Sham orbitals. Thus, the 4/ 
orbitals are localized as follows. Because of the inversion 
symmetry of the complexes, the 4/ orbital part of each 
Kohn-Sham orbital ipi is decomposed into the antisym¬ 
metric and symmetric parts: 

\iTfl = 4 (|1) + |2)) Ca,^ + 4 (|1) - |2)) (S5) 


where, |1) and |2) indicate the if xyz orbitals on the first 
and the second lanthanide sites, respectively. The abso¬ 
lute values of Ca,i and Cs,i for the occupied Kohn-Sham 
orbitals for the up spin part are shown in Fig. SI. As 
the antisymmetric and the symmetric levels, we averaged 
the Kohn-Sham levels: 


e/.a 


E occ. f-,2 
i 


E occ. 

i 


Ci 




E occ. 
i 


n2 


E OCC. ^2 

i ^ s,i 


(S6) 


In Eq. (S6), the sum is taken over occupied Kohn-Sham 
orbitals. With the use of the levels, the parameters t and 
A are derived (Table I in the main text). The transfer 
parameter is gradually decreasing as the increase of the 
atomic number because the ionic radius of the lanthanide 
shrinks. 


B. Calculation of it* —> 4/ electron promotion 
energy for 1 


The high- and low-spin states of the complex 1 were 
analyzed based on the Hubbard Hamiltonian: 


H = 

^ ^ ^ ^ ^ ^ ^7r*'^7T* 

'■ G 



2=1 ,2 -ycr a 



+ 

(7 


•o-C 27 (t^ 

+ 

^ ^ ^ ^ a' 

71 tt*'I'^' 77 * 4 ' 



2 = 1,2 ( 7 ( 7 , 7 'cr') 



+ 

^ ^ ^ ^ ^ ^ '^^ 7 cr'^ 7 r*cr^; 


(S7) 


2=1,2 7cr < 7 ' 



where 7 is the component of the 

4/ orbital, Uf 

and Utt* 


are the intrasite Coulomb repulsions on Gd and N2 sites, 
respectively, and v is the intersite Coulomb repulsion be¬ 
tween the Gd and N2 sites. 

The high-spin state with the maximal projection is de¬ 
scribed by one electron configuration: 

|lt,7r*t,2t), (S8) 


where the subscript “a” and “s” indicate antisymmetric where 1 and 2 are the lanthanide sites and f and j, are 
and symmetric orbitals, respectively. Comparing these spin projections. The if electrons which are not in the 
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FIG. SI. Contributions of the antisymmetric \Ca,i\ (red) and the symmetric (blue) combinations of the '^fxyz orbitals to 
each Kohn-Sham orbitals for the (a) Gd, (b) Tb, (c) Dy, (d) Ho, (e) Er complexes. 


4/a;yz orbital are not explicitly written here. The total 
energy ilps is 

= Eq + {2n + l)e/ + A + 2nv + n{n — l)u/, (S9) 

where Eq is the total electronic energy except for the 
electrons in the 4/ orbitals and tt* orbitals, and n is the 
number of the 4/ electrons in Gd^^ ion. For the low-spin 
state (tii, t type), the basis set is 

{|lt,U,2t),|lt,^U,2t),|lt,2i,2t)}. (SIO) 

Here, the configurations with the electron transfer from 
the 4/ to the tt* are not included because these configu¬ 
rations do not contribute much to the low-energy states 
due to the large energy gap A between the 4/ and the 
TT* levels. The lowest energy is 

= Eq + {2n -I- l)e/ -I- n{n — l)uf 
+ - ^A + 2nv + nuf — ^(A -|- 2nv — nuf)^ + . 

(Sll) 

The energy difference between the low- and high-spin 
states are 


Ai? = i?LS — ^'hs 
1 ^ 


(S12) 


nuf — (A -b 2nv) — •^(A -|- 2nv — nuf)^ + 

(S13) 

(u - VU^ + , (S14) 


1 
2 


where 

C/ = nuf — A — 2nv (S15) 

is the (averaged) electron promotion energy. Eq. (S15) 
shows that (i) the energy gap A significantly reduces the 
promotion energy and (ii) the promotion energy increases 
with the number of the 4/ electrons n. Using the transfer 
parameter t derived from the Kohn-Sham orbital, energy 
gaps between the high-spin state and low-spin state, and 
Eq. (S14), the averaged promotion energy U is derived. 


II. AB INITIO CALCULATIONS 
A. Fragment calculations for Ln^^ centers in 1-5 

To obtain the local electronic properties of the mag¬ 
netic ions, ab initio quantum chemistry calculations 
(CASSCF/SO-RASSI) were performed using Molcas 
[SI]. In the calculations, one of the metal ions in the com¬ 
plex was replaced by diamagnetic lanthanum ion (La^“'") 
and the ligands for the La ion were reduced (Fig. S2). 
Two point charges (—0.5 e) were put on each N atom 
creating the N2 bridge, where e (> 0) is the elementary 
charge. The latter is to include the electrostatic poten¬ 
tial from the unpaired electron of bridge. The co¬ 
valent effect is included later (IL'f in the main text). In 
the CASSCF calculations, all 4/ orbitals of the magnetic 
site are included in the active orbitals. The spin-orbit 
coupling is included in the SO-RASSI calculation. In 
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FIG. S2. The LnLaNj” fragment used in ab initio calcu¬ 
lations. Hydrogen atoms are omitted for clarity. The right 
lanthanide ion was replaced by La in the ab initio calculations. 


TABLE SI. Contractions of the employed ANO-RCC basis 
sets for the ab initio calculations. 


Ln 

7s6p4d2flg 

Si 

4s3p 

Ln 

7s6p4d2f 

0 

3s2p 

N (N 2 bridge) 

3s2pld 

C 

3s2p 

N (the others) 

3s2p 

H 

2s 


the SO-RASSI calculations the following CASSCF states 
were mixed by spin-orbit coupling: for Gd, 1 octet, 48 
sextet, 120 quartet and 113 doublet states, for Tb, 7 
septet, 140 quintet, 113 triplet and 123 singlet states, for 
Dy, 21 sextet, 128 quartet and 130 doublet states, for Ho, 
35 quintet, 210 triplet and 196 singlet states, and for Er, 
35 quartet and 112 doublets states. As the basis set for 
the calculations, ANO-RCC was used. The contraction of 
the basis set is shown in Table SI. The Cholesky decom¬ 
position threshold was set to 5 x 10“® Hartree. The ob¬ 
tained SO-RASSI wave functions were transformed into 
pseudo spin states (or pseudo J states) [S2-S5] to an¬ 
alyze the magnetic data using SINGLE_ANISO module 
[S6]. 

The obtained crystal-field (CF) levels are shown in Ta¬ 
ble S2. In all cases, the lowest spin-orbit states are dou¬ 
bly degenerate (Kramers doublet for Ln = Gd, Dy, Er) 
or quasidegenerate (Ising doublet for Ln = Tb, Ho). The 
ground CF states |^) are decomposed into the sum of the 
ground pseudo J multiplets \JM) [S4, S5]: 

.7 

IV')= 5] Cm\JM). (S16) 

M=-J 

The coefficients Cm are shown in Table S3. The con¬ 
tributions of the multiplets with the largest projection 
(|M| = J) to the ground CF states are 94.2 %, 96.4 %, 
97.1 %, 91.7 %, 78.6 %, for Gd, Tb, Dy, Ho, and Er, 
respectively. For each ground doublets, the g-tensors are 
calculated (Table S4). The Er complex is not magnet¬ 
ically anisotropic as much as the other complexes (Tb, 


TABLE S2. The lowest spin-orbit levels of Ln centers ob¬ 
tained by ab initio fragment calculations (cm“^). The cova¬ 
lency effect is not included. 


Gd 

Tb 

Dy 

Ho 

Er 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.099 

0.000 

0.982 

0.000 

0.329 

141.153 

179.143 

87.999 

74.691 

0.329 

142.222 

179.143 

88.534 

74.691 

0.631 

288.590 

320.747 

130.818 

118.034 

0.631 

295.694 

320.747 

147.454 

118.034 

1.108 

401.446 

406.717 

167.052 

166.279 

1.108 

435.925 

406.717 

202.496 

166.279 


490.539 

470.573 

224.559 

212.344 


531.201 

470.573 

241.625 

212.344 


547.372 

531.942 

246.712 

262.700 


730.715 

531.942 

284.704 

262.700 


731.087 

623.187 

296.344 

295.345 



623.187 

323.988 

295.345 



749.919 

327.012 

396.290 



749.919 

385.730 

396.290 




386.659 



TABLE S3. \JM) structure of ground GF doublet on 
center in 1-5 

Ln^+ 

Gd 

M \Cm\ M 

Tb 

\Cm\ 

Dy 

M \Cm\ 

M 

Ho 

\Cm\ 

Er 

M \Cm\ 

-7/2 0.971 -6 

0.694 

-15/2 0.986 

-8 

0.677 

-15/2 0.887 

-5/2 0.001 -5 

0.005 

-13/2 0.019 

-7 

0.005 

-13/2 0.112 

-3/2 0.225 -4 

0.123 

-11/2 0.164 

-6 

0.162 

-11/2 0.321 

-1/2 0.004 -3 

0.014 

-9/2 

0.027 

-5 

0.044 

-9/2 

0.168 

1/2 

0.077 -2 

0.024 

-7/2 

0.023 

-4 

0.084 

-7/2 

0.214 

3/2 

0.002 -1 

0.008 

-5/2 

0.007 

-3 

0.056 

-5/2 

0.103 

5/2 

0.038 0 

0.009 

-3/2 

0.010 

-2 

0.039 

-3/2 

0.105 

7/2 

0.000 1 

0.008 

-1/2 

0.004 

-1 

0.033 

-1/2 

0.024 


2 

0.024 

1/2 

0.002 

0 

0.027 

1/2 

0.031 


3 

0.014 

3/2 

0.001 

1 

0.033 

3/2 

0.021 


4 

0.123 

5/2 

0.001 

2 

0.039 

5/2 

0.011 


5 

0.005 

7/2 

0.000 

3 

0.056 

7/2 

0.012 


6 

0.694 

9/2 

0.000 

4 

0.084 

9/2 

0.016 




11/2 

0.000 

5 

0.044 

11/2 

0.002 




13/2 

0.000 

6 

0.162 

13/2 

0.005 




15/2 

0.000 

7 

0.005 

15/2 

0.000 






8 

0.677 




Dy, Ho). This is because the multiplets \JM) with small 
M {\M\ < J) are mixed more than the other systems. 
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TABLE S4. The g tensors for the lowest doublets of Ln cen¬ 
ters obtained from the fragment calculations. The transverse 
p-factors for Tb and Ho are zero because of the Griffith’s the¬ 
orem [S7]. 



Gd 

Tb 

Dy 

Ho 

Er 

gx 

0.492 

0.000 

0.0026 

0.000 

0.163 

9 Y 

0.824 

0.000 

0.0040 

0.000 

0.227 

gz 

13.439 

17.675 

19.6459 

19.422 

16.528 


B. Calculation of atomic J-multiplets of ions 

The excitation energies of the intermediate virtual elec¬ 
tron transferred states were replaced by the excitation en¬ 
ergies for isolated Ln^+ ion (Ln = Gd, Tb, Dy, Ho, Er). 
To obtain the energies, the CASSCF and the SO-RASSI 
calculations were performed with ANO-RCC QZP basis 
set [SI]. As in the case of the fragment calculations, all 4/ 
orbitals are treated as the active orbitals of the CASSCF 
calculations. In the SO-RASSI calculations, the follow¬ 
ing LS terms are included: for Gd^’*', ®P, ®F, for 

Tb2+, ^F, ^G, for Dy2+, ^F, ^G, for , 

and ^F, for Er^+. The excitation energies AF are 
shown in Table S5. 


III. ANALYSIS OF FIRST RANK EXCHANGE 
PARAMETERS 

As shown in Table II in the main text, the first rank 
part {k = k' = 1) oi the exchange interaction is isotropic 
Heisenberg type in all complexes, i.e., 

f7i±ii=Fi = —Jww 7^ 0; (S17) 

and the other J\q\qi are zero. The reason can be un¬ 
derstood analyzing the formula of the exchange interac¬ 
tion. The exchange parameter between J multiplet and 
isotropic spin 1/2 (Eqs. (2), (3) in the main text) is 


written as [S8] 


rr ^ ^}kqk'q' QajJk'xk^k' ia^o\ 

Jkqk'q' — 2 _^ 2 ^ — , A rnn-l-1 ’ (blOj 


X ajJ 


Uo + AF"+i 


where 


mm' ^ 

^ ^m-K'^^-K’^m'^l^m'kq^k' — q'lim’’ (S19) 


trmr* F the electron transfer between the 4/ with compo¬ 
nent m of orbital angular momentum and the tt* orbital 
of N2, = 3 is the magnitude of the atomic orbital angu¬ 

lar momentum for / orbital, x {li — k' < x < li + k') indi¬ 
cates a rank, ^ = -x, -a;-bl, ...,x, CfEfcg and Glf_^,i^^ 
are Clebsch-Gordan coefficients [S9], aj and J are the 
LS-teim and the total angular momentum of Ln^+, re¬ 
spectively, is the excitation multiplet energies of 

Ln^+, and G^^'jjk'xk and are functions of their sub¬ 
scripts. For the detailed description of x, G]2ljk'xk^ and 
see Ref. S8. 

Since the dependence of the exchange parameter (S18) 
on q and q' appears only in {t x (S19), the condi¬ 

tion for the isotropy of Jkq\q' is revealed from the equa¬ 
tion. First, we consider the cases where only the transfer 
between f±mo orbitals (mg = 0,1, 2,3) and the isotropic 
spin is nonzero for simplicity. The values of {t x 
are tabulated in Table S9. We find that the condition 
(S17) is fulfilled when mg = 2, while it is not for other 
mg. In the case of mg = 1, the nonzero terms with 
q = q' = ±1 are also the source of the anisotropic ex¬ 
change. When more than one set of / orbitals mg con¬ 
tribute to the electron transfer, the exchange interaction 
becomes always anisotropic. Finally, since Eq. (SI9) is 
independent of ions, the condition given above applies to 
the exchange interaction between any / electron ions and 
spin 1/2. 
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TABLE S5. Excitation energies with respect to the lowest J multiplet of isolated Ln^"*" ions (meV). 



Gd 



Tb 



Dy 


LS term 

J 

AE 

LS term 

J 

AE 

LS term 

J 

AE 


6 

0.000 


15/2 

0.000 

®7 

8 

0.000 


5 

182.104 


13/2 

307.374 


7 

458.212 


4 

333.856 


11/2 

573.765 


6 

859.148 


3 

455.259 


9/2 

799.172 


5 

1202.808 


2 

546.311 


7/2 

983.596 


4 

1489.190 


1 

607.012 


5/2 

1127.038 


6 

3412.346 


0 

637.362 


11/2 

1050.937 


5 

3756.005 





9/2 

1276.345 


4 

4042.388 





7/2 

1460.769 


5 

2369.357 





5/2 

1604.210 


4 

2655.740 




6p 

7/2 

4357.932 


4 

5667.150 





5/2 

4501.373 





Ho 



Er 





LS term 

J 

AE 

LS term 

J 

AE 





15/2 

0.000 


6 

0.000 





13/2 

638.260 


5 

850.945 





11/2 

1191.419 


4 

1197.911 





9/2 

1659.477 

^F 

4 

1560.065 





11/2 

3521.812 








9/2 

3989.870 







*F 

9/2 

2461.643 








TABLE S6. Kinetic contributions to the CF parameters Jkqoo 
(cm~^) for complexes 1-5. 


k 

q 

Gd 

kJ bgOO 

Tb 

Dy 

Ho 

Er 

0 

0 

-94.88 

-95.77 

-70.78 

-55.38 

-24.20 

4 

0 

5.86 xlO"® 

-30.11 

23.76 

10.00 

-11.15 

4 

±4 

3.50 xlO"® 

-18.00 

14.20 

5.97 

-6.66 

6 

0 

4.12 xlO"'^ 

-4.95 X 10"^ 

6.13 

-8.53 

4.29 

6 

±4 

-7.70 X 10"'^ 

9.27 X 10"^ 

-11.46 

15.96 

-8.02 


between J multiplets,” Phys. Rev. B 91, 174438 1-18 
(2015). 

[S9] D. A. Varshalovich, A. N. Moskalev, and V. K. Kher¬ 
sonskii, Quantum Theory of Angular Momentum (World 
Scientihc, Singapore, 1988). 
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TABLE S7. Energies of the CF multiplets (cm~^) of Ln cen¬ 
ters in 1-5 originating from the ground atomic J-multiplet of 
the corresponding Ln®^ ions, calculated with included kinetic 
contribution. Due to the latter, the ground CF multiplets of 
Ln centers are stabilized by 95, 110, 58, 55, 14 cm“^ for Gd, 
Tb, Dy, Ho, Er, respectively. 


Gd 

Tb 

Dy 

Ho 

Er 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.055 

0.000 

1.297 

0.000 

0.329 

168.191 

163.450 

95.958 

68.276 

0.329 

168.965 

163.450 

96.064 

68.276 

0.630 

316.136 

300.460 

129.220 

112.297 

0.630 

318.418 

300.460 

147.698 

112.297 

1.108 

426.550 

396.451 

174.642 

150.837 

1.108 

444.805 

396.451 

211.117 

150.837 


501.172 

458.147 

223.439 

197.730 


541.183 

458.147 

248.194 

197.730 


556.192 

510.902 

250.839 

249.831 


740.636 

510.902 

280.160 

249.831 


741.002 

604.121 

290.136 

276.581 



604.121 

320.445 

276.581 



747.503 

322.767 

387.217 



747.503 

371.677 

387.217 




371.948 



TABLE S8. Energy of the low-lying exchange KDs (cm in 
1-5 


Gd 

Tb 

Dy 

Ho 

Er 

0.000 

0.000 

0.000 

0.000 

0.000 

0.381 

207.619 

120.686 

105.154 

27.999 

0.643 

207.670 

120.686 

106.791 

28.000 

0.865 

210.623 

158.882 

108.718 

53.467 

1.187 

227.323 

164.275 

110.590 

64.209 

1.605 

362.446 

252.462 

146.181 

68.710 

2.112 

366.170 

273.941 

153.019 

86.241 

27.527 

369.751 

273.943 

160.514 

86.254 

27.761 

369.876 

293.589 

161.003 

99.987 
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TABLE S9. {t X for mo = 0,1, 2, 3. 


X 

gi 

52 


mo 






0 

1 

2 

3 

2 

0 

0 

3 1^12 |2 

7 I^Ott* I 

16|+12 |2 
\ 

10|+12 |2 
^l^27r* 1 

0 

2 

±1 

Tl 



10|+12 |2 

5 Ul2 |2 
“7p3ir* 1 

2 

±1 

±1 

0 

2/,12 \2 

7 V^zhl-TT* ) 

0 

0 

3 

0 

0 

0 


2| + 12 |2 
“3l^27r* 1 


3 

±1 

=Fl 


111+12 |2 

12 I^Itt* I 

2|+12 |2 
3l^27r* 1 


3 

±1 

±1 

0 

-WA.-f 

0 

0 

4 

0 

0 


15|+12 |2 

1 

6|+12 |2 
7F27r* 1 


A 

±1 

Tl 

5 1+12 |2 

3|+12 |2 

6| + 12 |2 

29ul2 |2 


"lil^OTT* 1 

~4l^l7r* 1 

~7p27r* 1 

28 1^371* 1 

4 

±1 

±1 

0 

3 /+12 \2 

“lil^ilTT* J 

0 

0 



